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Abstract 



The vortex velocity probability distribution for two distinct vortices is de- 
termined for the case of phase-ordering kinetics in systems with point defects. 
The n-vector model driven by time-dependent Ginzburg-Landau dynamics 
for a nonconserved order parameter is considered. The description includes 
the effects of other vortices and order parameter fluctuations. At short dis- 
tances the most probable configuration is that a vortex-antivortex pair have 
only a nonzero relative velocity which is inversely proportional to the distance 
between them. The coefficient of proportionality is determined explicitly. 
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I. INTRODUCTION 



It seems plausible that much of the structure one sees in the phase ordering of many 
materialsB! can be associated with the evolution and correlation among defects! like vortices, 
monopoles, disclinations, etc. These topologically robust objects grow out of the frustration 
suffered by a system with a continuous symmetry which is thermodynamically driven to 
align in a broken symmetry state. In the case of the n-vector model with the number of 
components (n) of the order parameter equal to the spatial dimensionality (d) one has point 
defects which are vortices for n = 2 and monopoles for n = 3. Because of the conservation of 
topological charge, the ordering in these systems is through the charge conserving process of 
vortex-antivortex annihilation. Topological constraints render the ordering in such systems 
to be largely independent of the microscopic details of the material. In this paper the 
following question is addressed: What is the probability, given a vortex at position f x with 
velocity v\, that one will find a vortex at position r 2 with velocity v-p- Clearly in answering 
this question, we obtain a tremendous amount of information about the dynamics of vortices. 

The calculation of the two vortex velocity probability distribution is a very involved 
process. In principle one could probe vortex dynamics by applying a force. Unfortunately in 
neutral systems it is very difficult to couple directly to the vortices. The two vortex velocity 
probability distribution serves this purpose by looking at the motion of one vortex in the 
fixed presence of another vortex a known distance away. 

The physical results of this calculation, carried out in detail for n = d = 2, are relatively 
simple to state. The appropriate probability distribution is a function only of the scaled 
velocities Ui = Vi/v for i = 1 or 2, and the scaled separation x — (fi — r 2 )/L(t). Here L(t) 
is the characteristic length in the problem which grows with time t after a quench as t 1//2 in 
the present case and drives the scaling behaviorS found in the problem. The characteristic 
velocity v, defined carefully below, is inversely proportional to L(t). For a given scaled 
separation x between two chosen vortices, the most probable configuration corresponds, as 
expected, to a state with zero total momentum and a nonzero relative momentum only along 
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the axis connecting the vortices: 

V\ = —v 2 = vx . (1) 

Moreover there is a definite nonzero value for v = v max for a given value of x. These most 
probable values are given as a function of x in Fig.l. The most striking feature of these 
results is that for small x the most probable velocity goes as 

Umax \^ ) 

where R is the unsealed separation between the vortices and k = 2.19 in dimensionless 
units defined below. The result giving v max inversely proportional to R is consistent with 
overdamped dynamics where the relative velocity of the two vortices is proportional to the 
force which in turn is the derivative of a potential which is logarithmic in the separation 
distance. Thus these most probable results are consistent with the short- distance behavior 
being dominated by the annihilation of vortex- antivortex pairs. From previous worki we 
know that there is low probability of finding like signed vortices at short distances. Thus 
our results giving the velocity as a function of separation should be interpreted in terms of 
annihilating vortex-antivortex pairs. The results for same signed vortices can also be carried 
out but is considerably more involved as discussed below. 

The work here builds on the work in Ref. [5] where the single vortex velocity distribution 
was determined. As in the single vortex case there are significant widths associated with 
these most probable results. The widths come about because of the existance of other 
vortices as well as fluctuations in the order parameter field. There are, as shown in Fig. 2 , 
significant widths in the probabilities due to the presence of other vortices and fluctuations 
in the order parameter field. For example, as x — > 0, while the most probable relative 
velocity is 4.38/i?, the half-width at half maximum for this quantity, in these same units, 
is 2.08/ R. In the large separation limit the probabilities become, as expected, uncorrelated 
and each has the distribution of velocities found previously^ for a single vortex. 

The analysis here is built upon previous work on the ordering kinetic of 0(n) symmetric 
systems. The best available theories! for the order parameter correlation function were built 



up in the early 1990 's and have led to the belief that we have a fairly good understanding of 
how to calculate the associated scaling function. It also has become clear that the order pa- 
rameter correlation function or structure factor is a rather structureless quantity which does 
not give a great deal of direct information about the underlying disordering agents. This led 
Liu and Mazenkoi to look directly at the correlations between defects in the scaling regime. 
The key new element in this work, as discussed in some detail below, was the realization 
that the positions of the vortices could be labelled by the zeros of the order parameter field 
which could, in turn, be mapped onto the zeros of an auxiliary field rh(x, t). They were able 
to show, following work by HalperinB, how one could write explicit expressions for the signed 
and unsigned vortex densities in terms of the auxiliary field fh(r,t). This then avoids the 
technically defeating step normally encountered which requires one to indentify the vortex 
positions. 

The signed vortex density correlation function was determined analytically in Ref( ^) 
in terms of the variance of the auxiliary field under the assumption the auxiliary field is 
gaussian. This calculation left the auxiliary field correlation function f(x) undetermined. 
Liu and Mazenko assumed that one could use f(x) determined from a treatment of the order 
parameter dynamics away from the defect cores. 

The charged or signed vortex autocorrelation function does not separate out all of the 
desired information since it mixes the correlations between like and unliked signed vortices. 
It is not difficult to introduce an unsigned vortex autocorrelation function. Between the 
signed and unsigned autocorrelation functions one can construct linear combinations which 
give the vortex-vortex and vortex-antivortex correlation functions. Unfortunately, for tech- 
nical reasons it is more difficult to determined the uncharged autocorrelation function. Only 
recently have these difficulties been overcome by Mazenko and Wickham@. They found the 
results, expected on physical grounds, that there is a depletion zone at short distances for the 
vortex-vortex correlation function indicating repulsion. Simulationsi and experiments^ also 
show a depletion zone at short distances for like-signed defects. This is expected on physical 
grounds since like-signed defects repel one another. There is a clear discrepancy between 



theory and simulation results at short-scaled distances. The theory shows a monotonic be- 
havior as the separation distance goes to zero. The simulation, however, shows a maximum 
at short separation distances and then falls rapidly to zero. The depletion zone seen in this 
case in the simulations is harder to understand physically since the pair is attractive and 
headed toward annihilation. While the theory satisfies the sum rule implied by topological 
charge conservation, it does not appear that this general constraint is satisfied by the sim- 
ulations. It appears that the short distance behavior in the simulations is contaminated by 
the choice of a vortex core distance which is comparable to distances associated with the 
unphysical depletion zone. 

It seems clear that it would be desirable to supplement this information on the spatial 
correlation of vortices with information concerning vortex velocities. It was recently shown 
by the authorll that one could write down an explicit expression for the velocities associated 
with point defects in terms of the order parameter field. A key ingrediant in this develop- 
ment is the identification of a continuity equation satisfied by the signed or charged vortex 
density. This continuity equation gives a fundamental expression for conservation of topo- 
logical charge in the system. Using the gaussian closure assumption one can determine the 
single vortex velocity distribution P[ui]. The most interesting physical result is that there 
is a large velocity tail which was interpreted there as arising from the high velocities in the 
late stages of vortex anti-vortex annihilation. Bray0 has used scaling arguments to obtain 
the same large velocity tail. The existance of these large velocities will be supported by the 
calculation carried out here. 

One common and concerning element in the calculations of defect correlation functions 
and defect velocity distributions is the requirement that the auxiliary field scaled correlation 
function, f(x), be analytic as a function of x for short scaled distances. For example the 
fourth order gradient (—^tf( x ))\x=o enters naturally into the analysis of Pf^i]. The need for 
analyticity in x for f(x) is not naturally consistent with the simplest self-consistent analysis 
of f(x) following a treatment of the order parameter correlation function. Mazenko and 
Wickham0 recently showed that one can construct the theory so that f(x) is analytic in x 



for small x, but this was at the expense of making the properties of the order parameter 
correlation function worseEl. 

The tension between using the order parameter dynamics to determine f(x) and the 
requirement that f(x) be analytic in order to treat defect dynamics has been, to a degree, 
releaved by the very recent work of Mazenko and Wickhamlll. They used the newly proposed 
continuity equation for topological charge to derive the equation satisfied by the auxiliary 
field correlation function under the circumstances that the field is constrained to be near a 
defect core. As discussed briefly in Sect.III.c of this paper, they find the clean result that 
the auxiliary field correlation function determined in this manner satisfies a linear equation. 
This result is self-consistent with the assumption that the auxilary field is gaussian. The 
solution of the associated linear equation has the Ohta, Jasnow and Kawasak0 (OJK) form 

f{x) = e"^ 2 (3) 

which is clearly analytic in the small x regime. They argue in Ref. |13| that the use of the 
gaussian assumption in determining defect dyanamics, such as -P(^i), has a stronger funda- 
mental justification than in the case of the determination of the order parameter correlation 
function. In the calculation of the two-vortex velocity probability distribution presented 
here it is assumed that the order parameter field can be replaced by a gaussian field in those 
portions of space near a vortex core and the associted auxiliary field correlation function is 
of the OJK form. 

II. ORDER PARAMETER DYNAMICS 

The system studied here has a defect dynamics generated by the time-dependent 
Ginzburg-Landau (TDGL) model satisfied by a nonconserved n-component vector order 
parameter ip(r,t): 
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where T is a kinetic coefficient, F is a Ginzburg-Landau effective free energy assumed to be 
of the form 

F = J Ag(W) 2 + K(H)) (5) 

where c > and the potential V is assumed to be of the 0(n) symmetric degenerate double- 
well form. Since only these properties of V will be important in what follows we need not 
be more specific in our choice for ff is a thermal noise which is related to T by a 

fluctuation-dissipation theorem. We assume that the quench is from a high temperature 
(Tj > T c ) where the system is disordered to zero temperature where the noise can be set to 
zero (ff = 0). It is believed that our final results are independent of the exact nature of the 
initial state, provided it is a disordered state. 

It is well established that for late times following a quench from the disordered to the 
ordered phase the dynamics obey scaling and the system can be described in terms of a 
single growing length L(t), which is characteristic of the spacing between defects. In this 
scaling regime the order-parameter correlation function has a universal equal-time scaling 
form 

(7(12) = <#1) • if (2)) = $F{x) (6) 

where ipo is the magnitude if} = \if>\ of the order-parameter in the ordered phase. Here 
we use the short-hand notation where 1 denotes (r^ti). The scaled length x is defined as 
x = (r*! — r 2 )/L(t) where L(t) ~ t 1 / 2 for the non-conserved models considered here. 

In previous work on the order parameter scaling function it was important to make a 
mapping of the order parameter if) onto an auxiliary field fh with the key requirement that 
away from defect cores 

f = ipom (A) (7) 

for the lowest-energy defects having unit topological charge. Physically one expects that 
near the defect cores 
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ip = am + b(m) 2 rh + ■ • ■ (B) (8) 

for charge ±1 defects where a and b are constants. In the theory for the order parameter 
correlations it is property A which is important. In the theory of defect motion, as presented 
here, it is property B which is important. In this paper only property B enters into the 
analysis since we always work near the defect cores. To complete the definition of the model 
one must specify the form of the probability distribution for the auxiliary field m. The 
simplest choice is a gaussian probability distribution for in with 

(m„(l)nv(2)) = C (12). (9) 

The system is assumed to be statistically isotropic and homogeneous so Co(12) is invariant 
under interchange of its spatial indices. In the scaling regime at equal times (t± = t 2 = t) 
we introduce the auxiliary field autocorrelation function mentioned in the introduction 

f(x) = C (f 1 t,f 2 t)/S (t) (10) 

and So(t) = Co (11) grows as L 2 (t) with time after the quench. 

III. TOPOLOGICAL DEFECTS 
A. Densities 

It has been emphasized in Refs.( H) that the signed or charged point (n=d) defect density 
can be written in the form 

p(R,t)=6$(R,t))V{R,t) (11) 

where the Jacobian obtained with the change of variables from the set of vortex positions 
to the zeros of the field ip is defined by: 

V(R,t) = ^e W)M2j ... ) ^e Vlji/2i ... )1/ „V m V' i , 1 V M2 '0^....V Mri ^„ (12) 
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where e Ml)At2 ,... )Atn is the n-dimensional fully anti-symmetric tensor and summation over re- 
peated indices is implied. The key point is that the zeros of the order parameter ip locates 
the positions of the vortices. The unsigned density, n(R, i), is given by 

n(R, t) = 6(tf(R, t))\V(R, t) \ . (13) 
The charged vortex correlation function is given by 

C pp (K,t)=<p(K,t)p(0,t)> . (14) 

While the unsigned vortex correlation function is given by 

C nn (R,t) = (n(R,t)n(p,t)) . (15) 

It is shown in Ref . ( f|) that the vortex- vortex and vortex-antivortex correlation functions can 
be expressed in terms of C pp and C nn . C pp was evaluated in Ref.( ||) using the gaussian closure 
approximation. As shown in Ref.( f|), the evaluation of C nn in this same approximation is 
technically much more difficult than the calculation of C pp because of the absolute value 
sign in the definition of the unsigned defect density n. 



B. Conservation of Topological Charge 



It was shown in Ref.( ||) that the charged vortex density satisfies the continuity equation, 



dp - V 



(K) 



(16) 



where 



(n-1) 



" e a,At2 , ■ ■ ■ ,Mn , ■ ■ ■ ,V n K Vl V^ 2 ljj U2 . . . . V^ n 1p Un 



(17) 



A key point here is that is multiplied by the vortex locating 5-function. This means 
that one can replace K in J^ K > by the part of K which does not vanish as ip — > 0. Thus in 
the case of a nonconserved order parameter one can replace Jg in the continuity equation 
by 
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jf = ^^ C 1 ^, e/?,M2,.,^^i,^,.,^V 2 ^iV /J2 ^ 2 ----V Mn ^ n . (18) 



Because of the standard form of the continuity equation Eg .(Jig), it is clear that one can 
identify the vortex velocity field as 

J® 



V 



(19) 



This form for the velocity field is used inside expressions multiplied by the vortex locating 
^-function. 

C. Use of Topological Charge Conservation to Determine the Auxilliary Field 

Correlation Function 

In previous workffl a rather successful scheme has been developed for evaluating the 
order-parameter correlation function, T[x) and, in turn, the auxilliary field correlation func- 
tion f(x). As indicated in the introduction this leads to the problem that the auxiliary field 
correlation function is rendered nonanalytic as a function of x for small x. Mazenko and 
WickhamE! have recently shown that this problem can be addressed in a different way. 
Rather than using the order parameter equation of motion to determine order-parameter 
correlation function they used the continuity equation for topological charge to determine 
the auxiliary field correlation function. As in the rest of this paper we use the fact that in 
quantities like p and v, we can replace ip — > m everywhere. Then we can determine f(x) by 
satisfying 

^ < p(l)p(2) >= Vf x) < 5m))jf\l)p{2) > +Vf 2) < p(l)5($(2))jf (2) > . (20) 

under the assumption that m is a gaussian field. The calculation of the left-hand-side of 
Eq.(pO|) amounts to the evaluation of C pp . This calculation was carried out in Ref.( P) 
and is straightforward since C pp factorizes into a product of gaussian averages which can be 

(2) 

evaluated using standard methods. The calculation of the average over Jp can be organized 
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in a similar fashion. In the scaling regime, after an impressive set of cancellations, one finds 
the rather simple result that 

C(2) 

-^xf' = V 2 f + n / (21) 

a 

where 

= J_ < (Vm) 2 > (22) 
n z 

S = aL 2 . (23) 

and we introduce the constant 

"=i • (24 » 

This equation for / is linear and has the simple solution of the OJK form 

/ = e-t* 2 (25) 

with the conditions 

n = (-V 2 /) \ x=0 = nn . (26) 

For simplicity we set \i = 1 which amounts to choosing L(t) = 2\^Tct and results in the 
result for / given by Eq.fl3|). 



D. Vortex Velocities 

As an important application of the result Eq. (|19| ) for the vortex velocity field v consider 
the velocity probability distribution function defined by 

n P(£i) = (nSfr - v)) (27) 

where v± is a reference velocity, v is given by Eq.(JIP|), n is the unsigned defect density, and 
n o — ( n )- was found in Ref.( |^) to be given by 

11 



r(f + 1) 1 

P (^l) = (7 ^2W2 7 TKT2)72 ( 28 ) 



where the the characteristic velocity v is defined by 



v 2 = (rc) 2 |^ (29) 



where is given by Eq.(|22|) and 



S 4 = I < (V 2 m) 2 > Jj^pl . (30) 
Using the OJK form for f(x) we obtain = a, §4 = ^ and = 



IV. CALCULATION OF THE TWO VORTEX VELOCITY PROBABILITY 

DISTRIBUTION 

A. General Development 

The main quantity of interest in this paper is the two velocity correlation function defined 

by 

C nn P[v h v 2 ] =< n(l)<f(i?i - v{l))n{2)6{v 2 - v{2)) > . (31) 

where V\ and are external labels while the v{i), for i = 1, 2 , is expressed in terms of the 
order parameter field via Eq.(|l9|). C nn P[v\, v 2 ] is normalized such that the integrals over 
v\ and V2 gives the unsigned defect density correlation function C nn which was determined 
previously in Ref.( f|). 

The first step in the evaluation of P[vi,V2] is to notice that we can replace ip by fh in 
the expressions for the unsigned vortex density and the velocity. Next we need to show that 
it can be expressed in terms of an average over a reduced probability distribution. In the 
Appendix we introduce the fields 

n 

W&A = 5(m(i))6(b(i) - V 2 m(*)) ]J - Vjpm^i)) (32) 

fl,V = l 
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which have the normalization 



„ n 

jd n b{i) J] dC^W^b] = 6(m(i)) . 



(33) 



Using this result we can insert the factors of W±W 2 into the expression for P[v\,v 2 \ and 
use the properties of the ^-function to replace all gradients and Laplacians of rh with the 
associated values constrained by the multiplying 5-function to obtain 
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C nn P[v u v 2 ] 



where 



n 

i=i 



d n b(i) n dc»(i)\v(t(i))m-mi)Xi))) 



G 2 (£,b) =<W 1 [C,b]W 2 [C,b] > 



(34) 



/ i. 



(35) 



with 



4 2) (e(o,6(0) = 



r c 



n-l ! 



(36) 



(37) 



— * 

The gaussian average giving G 2 (£, b) is worked out explicitly in the Appendix. In the course 
of this calculation it is required that one make a change of variables from to a new set 
given by 



C(0 = 



(38) 



and i?^ is an othonormal matrix with the additional property that Rjp = where is 
the unit vector pointing from vortex 2 to vortex 1. Since det(R) = 1 the change of variables 
from £ to t is simple, 



n<(j) = n^(j) 



(39) 



I'M? 
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and 



mm = vim 



(40) 



The one place in this change of variables where one must show care is for the current J( 2 \ 
We have 



Tc 



(n-1)! 



^a,/f2,...,/in^fl,^2,---> I> 'ra 



Tc 



Clearly if we multiply this expression by R 1 ^ 1 and sum over a we obtain 



(41) 



Tc 



(n-1) 



— (det i?) e /3l , / 3 2 ,..., /9n e i , 1)I , 2i ... );/ „fo l/1 (i)^(i),...,^(i) 



4 2) (*(0,?(0) 



(42) 



Multiplying by i?^ 1 , summing over /3i and using the orthonormality of the matrix R gives 



(43) 



Because of the rotational invariance of the <i-dimensional 5-function we have 



5{v t - v(£(i), &(«))) = 6(u(i) - v(t(i), b(i))) 



(44) 



where 



(45) 



Thus the u — 1 component of is the longitudinal component along R. We then have that 



r 2 

c nn p[v 1 ,v 2 ]= j n 

i=l 



fl,V=l 



G 2 (t,b) . (46) 
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The next step in the analysis is to perform the integration over the b variables. Toward 
this end we use the representation 

6(u(i) - v(t(i),b(i))) = I ^^ e -i<i)<i)^mi)ki))-m ( 47 ) 

J (27r) n 

and we make explicit the b{i) dependence by writing 

v{t(i),b{z)) -z(i) = a v {i)b v {i) (48) 



where 



Tc 1 



z a (i)N va (i) (49) 



and 



r c i 



N ™{*) = - JZ ^^u^^ ^---^uu 2 ...uj%(i)-.r^ (50) 



Next one must make the b dependence of G2(t, b) explicit. We have from the Appendix that 

G 2 ^,b) = G T (t T )G L (b,t L ) . (51) 

The transverse part of G 2 does not depend on b(i), while the longitudinal contribution can 
be written as 

1 1 6 ~* ~* 

nit T \ _ _ t p-lEc^i^-MAf-'W (rn) 

l*L{b,t L ) (27r )3n( detM )n/2 e " W 

— * 

Where the matrix M is discussed in detail in the Appendix and the h a are defined by 
Eqs.( |A47| - |A51j ) . Using the explicit expressions for the h a we can write 



£ K ■ h^M-Xp = S b J2b(i) 2 + 2C b b(l) ■ 6(2) + 2^6(z) • S(i) + S° L + 2C° L t L (l) ■ f L (2) 

a, /3=1 i i i 

(53) 

where we have defined 
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C 6 = (M- 1 ) 34 = (M- 1 )^ 



5(1) = (M-^astKl) + (M- 1 ) 36 r L (2) 



5(2) = -(M-^aet^l) - (M" 1 ^^) 



66 



Cl = (M-^se = (ikT 1 ) 
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(54) 
(55) 
(56) 
(57) 
(58) 
(59) 



The matrix inverseM 1 is also discussed in detail in the Appendix. It is convenient to define 



G° L (b,t L ) 



1 



1 



(2tt) 3 ™ (det M)«/ 2 
and the 6-dependence is explict when we write 



(60) 



G 2 (e, b) = G r (t r )G°jr L ) e -M^E l ^) 2 +2c b 6(i). ; ;(2) + 2E l K i )-s«] 



(61) 



Putting this together we have 



c nn p[v 1 ,v 2 ]= j n 



n « 



l^=i (2tt)» 
where the 6(i) integrations are isolated in 



z>(*(0)| 



G r (t r )G'°(r i )e-^ i ««-^)j (62) 



J = [ TT<f»&(i) e *£i&(0-[a«^«] e -M*^ 
-> i 

This integration is of the standard gaussian form. If we define 

A{i) = a(i) + iS(i) 

then we have the result 



(63) 



(64) 
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where 



./„= [fj e~**i 



Q 



(65) 



Q = J2A 2 (t)-2f b A(l)-A(2) 



(66) 



and 



7, = (i-A 2 )- 1/2 



(67) 



fb — Cb/ Sb ■ 



(68) 



The next step is to do the z(i) integrations. We can highlight the z dependence if we 
remember that 



a v (i) = z a (i)N va (i) . 



(69) 



It is then a matter of straighforward algebra to show that 



Q = - E ^(0 + 2/^(1) • 5(2) + E E ^(i)E a ,(ij)z,(j) 

i ib a/3 ij 



where 



E al3 (ij) = -£-n a p(ij)[5ij - fbSj,i+i] 



and 



(70) 



(71) 



(72) 



Here we have introduced the convenient notation that the index i is periodic, so that if i = 2 
then % + 1 = 1. Putting these results together we have 

2 



C nn P\v x ,v 2 \ = ff[ 



II dt;(i)\v(t(i))\ 



G T (t T )Gl(t L ) 



x 



( 27r )» ^y e ^E l ^«-2/ 6 5(l).5(2)] Ji 



(73) 
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where the z integration is given explicitly by 



where 



i=l 



d n z(i) 
(2tt)™ 



(74) 



— * — » 

C/j = u(i) + d(i) 



(75) 



and 



J\ is again of the standard form for a gaussian integral so 



-.exp 



l -U a {i){E- l U{ij)Up{j) 



(76) 



(77) 



(2tt)«/2 

where, again, we need the determinant and the inverse of a matrix, in this case E. Let us 
look at the inverse first. If we note the important result (used in Ref.( |5|)) 



N aa {i)t» a {i) = -TcS 



(78) 



where we do not sum on i, then 



n a/3 {ij)t a Jj) = -TcN Q 



(79) 



where we do not sum on j. These identities suggest that we try a solution for E 1 of the 
form 



(80) 



with Cij to be determined. Inserting this ansatz into the equation defining the inverse we 
easily find that 

1 



(r c y 



(81) 



In computing del E we use the fact that 



det E 



det (E 



-v 



(82) 



and that E can be written as the matrix product 



(83) 



vv'kt 



so that 



det E 1 = det t u a {i)5ik det (eke) <W det Sijta (j) 
= V{t{l))V(t{2))(det e) n V(t(l))V(t(2)) 



(84) 



and 



det e 



(To 



1 1 S" 2 

(S 2 — ( '~ s - 



(rc) 4 7b 2 



(85) 



Pulling all of this together leads to the result 

1 



C nn P[vi,v 2 } 



(Tc) 



2f< 



2 n 

n n ^w^w)) 



G T (t T )G° L (t L ) 



xe 2 s b 



3^tEi^W-2/^(l)^(2)] 



1 



£ u^HE-^MUpV) 



ija/3 



(86) 



and the important point is that one does not have an absolute value sign left in the Jacobian 
factors. Turn next to the argument of the exponential in the last line of Eq . (|36"D . After a 
substantial amount of algebra we find 



J2S 2 (i)-2f b S(l)-S(2) 



+ £^(i)(ir%(ij)W) 

ija/3 



Ev«-v(j>ii- F Ev(i)'* 
ij 1 r i 



(87) 



where 



v„(0 = 



(88) 



Putting this together leads to 



1 



1 



\n(n-l) 



(27r) 3n (rc) 2 " (det M) n / 2 \2tiS t 
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J 2 



(89) 



where the final integration is over the matrices t^(i): 

2 n 



■h 



n n ^(o^ 2 (*(o) 



1 = 1 fl,V = l 



e 2 



\A(t) 



(90) 



where 



2 n n 



>T a=2 v=\ Li=l 



Y\m\ 2 - v<{\)ti{2) 



+2(^4(1) • 4(2) + £V(i) • V(j>^ - =-EV(») • S{i) 

ij i 



(91) 



A(t) is clearly a quadratic form in the matrix t v Ji). The quantities St, /t and 7t govern 
the transverse modes and are defined by Eqs.( [A56l - |A59|) in the Appendix. After sufficient 
rearrangement A(t) can be written in the final form 



(92) 



where the matrix W plays a central role in the theory and is given by the manifestly sym- 
metric form 



w aP {ij) = 6 aP d^ + u a (i)n$u P ti) + u a {i)nf)L p {j) + L a {i)d§ Up {j) 



(93) 



where 



d a = 5 a , L S L + (1 - S a>L 



± 
St 



(94) 
(95) 



d c a = 5 a , L C° L + (1 - S a , L )f T 



± 
St 



(96) 



(i) _ S b C b 

lj ~ v (r c y + i+1,J (r c ) 2 



(97) 



Tc 



Te- 



rn 
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and 



L a (i) = 5 ayL f]i 



(99) 



where 



vi = -m = 1 



(100) 



The final integration over the matrices t v a {i) is not of the standard form evaluated so far but 
instead there is the polynomial V 2 (l)V 2 (2) multiplying the gaussian in the integrand. It is 
technically important that there are no absolute value signs left in this expression and the 
integral can be evaluated by introducing a field g^i) which couples to t v a {i) via 



^(0->-^(*) + E£(0£(0 • 



101) 



If we consider 



Mg) 



2 n 

n n wawm) 

1=1 /U,I/=1 



e 3^(*) e E wi ^«tS(0 



(102) 



then any polynomial can be generated by taking derivatives with respect to g. Using the 
explicit expressions for the V 2 {t{i)) we have 

d d A d d 



n q q n d 8 



(103) 



where 



■Ms) 



2 n 

n n w 

1=1 li,v=\ 



e 2 



5^(*)pEa«^«*S(< 



(104) 



J 3 ((yf) is now of the standard form and we have 



M9) = WWy^ exp 



ija(3v 



(105) 



where A a p(ij) is the matrix inverse of W. It is then straightforward to take the derivatives 
with respect to g and then set g to zero to obtain 
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X 

<J=1 



J * = mw^- J > (106) 

where 

n n n n 

|Ul.../i n =l ^...^=1 v^...v n =\ v' x ...v' n =\ 

n 

II [A m ^(11)A^ ct (22) + A^, (12)A^„„(12) + A MctIV (12)A m , (12)] (107) 

cr=l 

and the final result is 

1 1 1 / 7 T \ «(«-!) (27r) n2 

C„„P[Vi, V 2 ] = (27r)3 „ (rc)2 „ M)n/2 {-^) {def W) n,2 J ± ■ ( 108 ) 

Finally all of the integrals have been evaluated. What is left is to evaluate the determinant 
and matrix inverse of W . 

B. Large x Limit 

As a check on the preceeding analysis it is useful to work out the large scaled distance 
limit where we expect the probability distribution to factorize into the product for each of 
the tagged vortices. In this limit, using the results from Table I, we find that the matrices 
entering W are in a diagonal form 

d% = D5 t3 (109) 

where 

D = - (110) 
a 

and 

ng> = o • (ii2) 
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The matrix W can then be written in the partially diagonal form 



W a p(ij) = Sij [5 a pD + u a (i)up(i)\ 



where 



LuJi) 



(113) 



(114) 



Yc\j2no 

Notice that there is no longer a difference between the longitudinal and transverse directions 
as expected. The inverse matrix A then satsifies the equation 



DA a/3 (ij) + u a (i) u^(i)X^(ij) = 5^5^ . 
This equation is in the form of a trap which can first be solved to obtain 



and the full inverse is given by 

A a /3(u) = Sij 



a/3 



U a {l)Up{l) 



D 2 



1 + Euul(t)/D 



(115) 



(116) 



(117) 



With these results it is easy to see that the quantity J4 can be written in the simple form 

J 4 = {nlfdet A(ll))det A(22) (118) 



and 



det W = det W(ll))det W(22) 



det A(m)) 



(119) 



(120) 



det W{ii)) 
It is then easy to see that 

det W[ii] = D n [l + v 2 Jv 2 } (121) 

where v 2 is given by Eq. (p9"D . If we then carefully keep track of all the factors we see that 
Eq. ( |108j ) reduces to 

\imC nn P[v u v 2 } = noPK]noP[tf 2 ] i 122 ) 

where n P[vi] is given by Eq.(30) in Ref.( 0) and, after proper normalization, leads to the 
expression for the single vortex velocity probability distribution given by Eq.(^). 
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C. n=d=2 Case 



The general expression for C nn P[vi,v 2 ], is complicated. Let us restrict ourselves here to 
the case of n = d = 2 where det W and J 4 can be evaluated explicitly. Let us define 

det A(ij) = A 11 (ij)A 22 (ij) - A 12 (ij)A 21 (ij) (123) 

where the matrix A a p(ij) is either W or its inverse A. We can also define 

Q A = L4 n (22)A 21 (21) - A 21 (22)A n (21)] [A 12 (11)A 22 (12) - A 22 (ll)4 12 (12)] 

- L4 12 (22)A 21 (21) - A 22 (22)A n (21)] [A 12 (11)A 21 (12) - A 22 (ll)^ n (12)] 

- L4 n (22)A 22 (21) - A 21 (22)A 12 (21)] [A 11 (ll)^ 22 (12) - A 21 (ll)4 12 (12)] 

+ [A 12 (22)A 22 (21) - ^22(22)^12(21)] [A n (H)A 21 (12) - A 21 (ll)A n (12)] . (124) 

In terms of these quantities we have 

det W = det W(U)det W(22) + det W(12)det W(21) + Q w (125) 

and 

J 4 = 4 [det A(ll)det A(22) + Met k{Y2)det A(21) + Q A ] . (126) 

It is clear that the last nontrivial step before evaluating C nn P[vi,v 2 ] is to determine A a/ g(ij). 
This will be carried out in general in section IV. E , however most of the important physics 
can be extracted in the problem by considering the simple case where the transverse velocities 
are both zero. In this case one can make substantial analytical progress. 

D. Zero Transverse Velocities 

Before tackling the complete determination of P[vi, vi] it is very instructive to study the 
much simpler case where the transverse velocities are set to zero. This case is of interest not 
only because it is simple but also because it is the most probable situation. The most likely 
situation is that each of the two tagged vortices will have zero transverse velocity. 
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If the transverse velocities are zero then the problem simplifies since the matrix W 
reduces to the diagonal form 



W aP (ij) = 5 a pD« (127) 

where 

Dfj = d% + 5 a , L [u(i)nW U (j) + umt^m + Vi^u(j)] (128) 

and u(i) = ul{i)- Clearly the longitudinal and transverse degrees of freedom are uncoupled 
and we have after some manipulations 

Dl = d T [Sy + f T 5 jti+1 ] (129) 

while 

Dfj = dibij + bS jji+1 (130) 

where 

a l = (M-% 5 + (M~% 3 u(t) 2 -2r ]t (M-% 5 u(i) , (131) 
b = (M-^se + (M- 1 ) 34 m(1)m(2) + (M-^ae [«(2) - u(l)] (132) 

and 

«(0 = ^ • (133) 

In this case we see that W is a relatively simple matrix. We need its determinant and then 
its inverse on the way to evaluating the quantity J4. The matrix of W for general n is simply 
given by 

det W = (det D T ) n ~\det D L ) (134) 

where 
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det D T = (d T ) 2 (l-f 2 ) 



Q2 n'2 



(135) 



and 



det D L = a x a 2 - b 2 



(136) 



The quantity det D L is key in the development and we shall return to it soon. First we 
need to evaluate the inverse of W to complete the calculation. In this case this involves the 
solution of the equation 



k 



which is easily found to be given by 



Aa/?(u) = 5 af3 [5 a , L A L (ij) + 5 a>T A T (ij)} 



(137) 



(138) 



where 



Al( ^ = det D L ^ ai+l6ij ~~ 



and 



(139) 



A r(ij) = def DT d T [8ij - f T Sj,i+i\ ■ 



(140) 



Using these results one can work out the quantity J 4 for the case n = d = 2 with the result 



J 4 = 4 



det A(ll)det A(22) + 3 (det A(22)) 2 - A L (22)A L (11)A|(12) - A T (22)A T (11)A^(12) 

(141) 



where we have used the fact that the matrices Ai(ij) and Ax(ij) are symmetric. Putting in 
the explicit forms for A we obtain 

4 



J A = 



det D L det D T + 2b 2 (d T f T ) 2 



(det D T ) 2 (det D L f 

Putting all of this together for n = d = 2 and the transverse velocities zero we have 



(142) 
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C nn P[vi,v 2 ] 



1 



(Tc) 6 det M 



1 



+ 



2b 2 C 2 



(143) 



(det D L ) 2 det D T (det D L f 

Clearly the next step is the explicit evaluation of det D L . Using the expressions for 
and b given by Eqs. (|131|) and (|132|) we obtain, after some algebra, that 

r 4 



det D L 



a 



D D E 



7o + ^(«(1)-«(2))+7a(« 2 (1)+« 2 (2)) 



73 ~, 



+^ B u(l)u(2) + |m(1)5(2) (u(l) - u(2)) + 74 n 2 (l)n 2 (2) 



where the scaled longitudinal velocities are defined by 



(144) 



u(i) = u(i)L 



(145) 



and 



(2) (2) 
70 = E 



(146) 



7i = 2 (l-/)4 1) 4 2) + (l + /)4 1) «: i 



o 



(147) 



ft f\ DE , (t , f\ D ° , ^ 2 
7A=(l-/)^ + (l + /)^ + T 



148} 



IB 



/)— + (1 + /)— 

4(7^ \ ^ ' 0-3 2 



72 = (1 - /)«g>«g } - (1 + /)*^ - 4(1 - f)K%K E 



(0)(2) 



1-2^(1)^(1) 



(149) 
(150) 



73 = 2(1 - f) 



K (i) K (o) _ (1) (0) 



(151) 



74 = (1 - / 2 )4 0) 4 0) 



(152) 



where the k's are given as functions of / by Eqs. (|A91j - [A96| ) in appendix A. Do and D E are 
given as functions of / by Eqs.( |A83| - |A84D in appendix A. Note the relationship 27a — 1b = 
72- Notice the crucial result that after rescaling the velocities by a factor of L, the time 
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dependence drops out of det D L . This will eventually lead to the result that scaling holds 
for the probability distribution at late times if we rescale velocities in this manner. 
We also need to express the quantity b in terms of the /t's. It is convenient to write 



b = b + b E 



where 



bo = 



AD 



o 



4 2) - (1 - /)^fi(l)fi(2) + (1 - fWo> (5(2) - fi(l)) 



~( )„~, 



(153) 



(154) 



a 



bE = ^ 



4 2) + (1 + f)^>u(l)u(2) + (1 + f)^> (fi(2) - «(!)) 



The last ingrediant needed to evaluate the probability distribution is 



Ck f{x) 



Ct = ~^ = -a 



X 



a f{x) 



(155) 



(156) 



using the OJK form for fix) in the last step. 

Let us look first at the small x limit. Since the OJK form for fix) is easily expanded 
in a power series in x and we can extract to leading order in x: 7o = 32x 4 , 71 = — 24x 5 , 



72 = fa: 6 , 73 = -fa; 7 , 74 = %1a = |^ 2 , 1b = 27a- We also need 



Dp* = 16x 2 a 3 



(157) 



D = <r A — ■ 



48 



(158) 



Notice that the 7a and 7_b dominate the expression for det D L in the small x limit and we 
can write to leading orders in x 



det D L 



7o + X4(m(1) + m(2)X 



D D E 

plus terms which are higher order in x. We can write this in the more convenient form 



(159) 



det D 



<? 7o 
D D E 



1 + v 2 /v* 



(160) 
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where after some algebra we obtain 



4- 

o_ 7o 



96 



(161) 



V = «(1) +u(2) 



(162) 



and 



\/ 2 = — = 8x 2 

1A 



(163) 



Similarly we find that for small x that b is dominated by bo and given to leading order in x 
hy b = — We also need Ct = u to leading order in x. The probability distribution is 
dominated in the small x limit by the term proportional to b 2 . The other term is down by 
a factor of x 4 . Putting all of this together we obtain 



a 



2\ e 



C nn P[vi,v 2 ] = — 



1 



v r c ; {i + v 2 /v 2 y 

Then as x — > we find, with increasing probability, that 

V = u(l)+u(2) = . 



(164) 



(165) 



This is just the physical statement that there is very low probability that there is a nonzero 
momentum of the center of mass of the two tagged vortices. Thus all of the action is in the 
center of mass where we can set 



u(2) = -u(l) = « . 



(166) 



If we return to the probability distribution for the case where the COM momentum is 
zero, then 

,4 



det D 1 



a 



D D E 



=r, „~,2 , zr, ~3 i - ~4 



70 + llU + 72« + 73« + 74« 



(167) 



and 



a 



^D 



4 2) + (1 - f)^>u 2 + 2(1 - f)K^U 



(!).-. 



(168) 
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r -4 2) - (1 + f)nfu 2 + 2(1 + . (169) 

In the small x limit these reduce to 

b = b + b E = -^f p (170) 



where 



u u 2 , _ . 

a = 1 -4 + « (171) 

96 

det{D ") = —f L (172) 

with 

= 1 — —it + — — w 2 — -w 3 + — — -m 4 (173) 
J 4 48 24 384 v ' 

and the scaled velocity is given by 

u = ux . (174) 
The probability distribution is given in the x — ► limit by: 



a 



2 \ 6 f2 



C nn P[v u v 2 ]= j^-j ^ . (175) 

A key conclusion we can draw at this point is that it is only the combination u = ux which 
enters the probability distribution with high probability as x — > 0. Thus the relative velocity 
increases as 1/x as x — > 0. We plot /|//| as a function of w in Fig. 2. We that the most 
probable values of the relative velocity as a function of x for small x are given by: 

_ Tc k Ten 

VL = Tx = ^r (176) 

with k = 2.19 a pure number. This is the result quoted in the introduction. 

If we then plot the two vortex velocity probability distribution function for zero transverse 
velocities in the COM for general x as shown in Figs. 3 and 4, we obtain the most probable 
relative velocities as a function of x as shown in Fig.l. The interpretation that this is the 
interaction between vortices and anti-vortices holds only out to modest values of x where 
the population of same signed vortices begins to appear (See Ref.( £|)). 
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E. General Evaluation 



The complete determination of the two- vortex velocity probability distribution as a gen- 
eral function of V\ and v-i can be carried out in the n = d = 2 case if one can invert the 
matrix W a p(ij) to obtain its inverse A a p(ij) defined by 

^2W a „(ik)A^(kj) = 8^ . (177) 

This inversion is a quite unpleasant task if one heads in the wrong direction. It is useful 
in order to make the development more transparent to introduce a mixed operator notation 
where W a p is an operator in the space associated with the indices i and j 

W afi (ij)=<i\W a p\j> . (178) 

Then the matrix W a p is given by 

W aP = 5 a/3 d a + uM {1) Uf3 + u a VL {2) L p + L a tt {2) u p (179) 

where u a and L a are diagonal operators 

< i\u a \j >= u a (i)5ij (180) 

< i\L a \j >= L Q (i)6ij . (181) 

The key idea is that if we can write W a p in the form 

W Q p = 5 aP D a + P a Pp , (182) 

where Pp is the transpose of P a , then we can carry out the inversion straight away. Let us 
first show this and then return to show that W can be written in the assumed form. 
We want to invert the equation 

E^W^'W • (183) 
Inserting the assumed form, Eq. (|182|) , for W we obtain 
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D a A al 3 + P a PfiA-ii/3 — S a p . 



;i84) 



Multiplying from the left by the matrix inverse of D a this becomes 



a/3 



(185) 



This equation is then in the form of a trap for the quantity Y,u P^nP ■ Multiplying Eq.( |185| ) 
by P a and summing over a we obtain a closed equation for J2n P^vP which has the solution 



186) 



where the 2x2 symmetric matrix Q is defined by 



q = J2p,d, 1 p, 



(187) 



This leads directly to the final result 

Kp = D-H^ - D- l P a [l + Qj^PpDp 1 (188) 

which is clearly symmetric. This gives a practical expression for the inverse once one has 
identified the matrices D and P. 

The key observation which allows one to write W in the desired form given by Eq. (|182|) 



is that the matrix Q)y can be factorized in the form 



where 



n 



(i) 



k 



(189) 



+i 



(190) 



and 



2Tc 



1 + h + Jl-fi 



(191) 



2Tc 



l + f b -Jl-fi 



(192) 
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Using this factorization result it is then easy to show that W can be written in the form 
Eq.(M with 

D ij = d % ~ <W Yl C L (ik)C L (kj) , (193) 

k 

P a {ij) = u a (i)Uij + C a {ij) (194) 

with 

Ca(v)=^a(0E n 2W • (195) 
k 

Combining these results one has an explicit expression for the two- vortex velocity probability 
distribution for arbitrary velocities. The major qualitative feature of including the transverse 
velocities is to allow one to look at the widths of the distributions in the transverse directions 
since we find the most probable configurations are those where the transverse velocites of 
both vortices are zero. These widths turn out to be comparable to those associated with 
the longitudinal modes. 



V. DISCUSSION 

In the analysis here we have looked at the correlation between vortices regardless of 
their signs. At short relative distances, where it is unlikely to have two vortices of the same 
sign, one can interpret the results in terms of vortex-antivortex dynamics. It is clear that 
one can go further, as discussed by Mazenko and Wickham@, and separate the probability 
distribution into that for vortex-vortex and vortex-antivortex pairs. The key idea , which 
is essentially equivalent to that used in the case of spatial correlations, that a factor of 
V+(l) = |(1 + sgn T>(1)) restricts one the positive charge vortex sector, while "P_(l) = 
|(1 — sgn P(l)) restricts one to the negative charge anti-vortex sector. Thus the probability 
for vortex-vortex correlations is 

C VV P VV (12) =< n(l)5(^ - v(l))V+(l)n(2)5(v 2 - v(2))V+(2)) . (196) 

33 



The vortex-antivortex contribution is given by 



C av P av {\2) =< n(l)5(v x - v(l))V-(l)n(2)5(v 2 - v(2))V+(2)) 



(197) 



These quantities can be multiplied out and, and after using symmetry to show that the 
correlation between the signed and unsigned quantities are zero, can be expressed in terms 
of the probability distribution determined in this paper and 



which has not yet been computed. It is expected that P pp (12) will be difficult to determine 
because of the addition factors of the sgn D. The analysis will be essentially identicle in 
structure up to eq.(125) with the expression for J 2 showing the replacement 



The resulting integral for J 2 can not then be represented in the product form given by 
Eq. ( |103j ) . This remains a problem to be solved. 

In principle Eq. (|108|) givens an expression which can be integrated over all velocities to 
give C u and determine the overall normalization. It is not clear how to do this analyti- 
cally since the velocities appear in a complicated fashion in det W and J4. A numerical 
determination is quite feasible. 

In this paper we have shown how one can make progress in an analysis of the dynamics 
of point vortices in the context of phase ordering kinetics. The results include the effects of 
other vortices and order parameter fluctuations on the dynamics of the tagged vortices. The 
results appear completely physical and the determination of the relative velocity as short 
distances appears to be a useful result. The method used here appears to generalize easily 
to the case of string defects. This will be the subject of subsequent work. 



C PP P PP {12) =< p{l)8{v x - v(l))p(2)5(v 2 - v{2))) 



(198) 



2> 2 (t(*))-P(t(z))|X?(t(z))| 



(199) 



APPENDIX: GAUSSIAN AVERAGE 



In this appendix we work out the gaussian average 
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G 2 (Z,b)=<W 1 [Z,b}W 2 [Z,b]> 



(Al) 



where the W's are defined by 



\x,v=\ 



(A2) 



where we have already assumed that n = d in the product. 

The first step in the evaluation of G 2 is to use the Fourier representation for the 5-function 
to obtain 



Wi[£,b]= / dQ\i]e l ^ m(t) e 



,;l R tf i .fR(i) R -iii.(6(i)-V?nl(i)) e -iE^ ^(0(^(i)-V M m„(i)) 



(A3) 



where we have defined 



dm = 7^7^ n 



dk;(i) 



(27T)-(27T)"^i 1 L 27T J 

We can rewrite this in the more useful form 



)m a (l) 



where 



/*-/ 



(A4) 



(A5) 



(A6) 



(5(lt) 



(A7) 



and V/*-, is the /Z^ component of the gradient acting on r*j. The average of interest can then 



be written as 



G 2 (f,&) = y dfJ[l]dfJ[2]< 



Eli s-5(0+£u„=i 



2 n 
Li=l a=l 17 



> 



(A8) 



The average is of the standard form for a gaussian average with the result 



< exp 



2 n 



i i J 



.i=l a=l 



(A9) 
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where 



2 n 



A) = -E E [ dl [ d2H^(l)Hf (2)00(1,2)5^ 

i,j=l a,/3=l 



and we have used 



(A10) 



< m a (l)m p (2) >=C (1,2)<^ 



Inserting the expression for H into A we need the following definitions: 



C (n) = S 



(AH) 



(A12) 



Vj)Co(y) 



5 (2) 



(A13) 



ll=J 



(A14) 



I • • = 5 S {2) 

\l=J u jJ,V>~> 



(A15) 



Using the fact that C (12) depends only on the magnitude of R — f\ — f 2 we convert all 

— * 

derivatives to those with respect to R: 



V{i } C7o(12) = C'oR, 



(A16) 



Vf 2) C (12) = -C' R, 



(A17) 



where the prime indicates a derivative with respect to R. Going further, for i — 1 and 2, 

(A18) 



V^C (12) = V 2 R C (R) = C> + — — , 



^fl)^(2)Co(l2) = —CqR^R u — -J£ (finis — -R M -Riy 



(A19) 



V? 2) Vf 2) C (12) = -p^ 



(A20) 
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Vf 1) Vf 1) C7 (12)=p4 (A21) 

where 

P = C^ + ^^-(CS-^) = (Vj t C (R)Y ■ (A22) 

We see that it is then natural to use the coordinate system parallel and orthogonal to 
R. Indeed we can introduce the orthonormal set Rp where 

± KK = V (A23) 

a=l 
n 

Y, = ■ (A24) 

The only other thing we need to know about this set is that 

Ri = Rp ■ (A25) 

Next we define 

^(0 = 1:^(0 (A26) 
n=i 

which can be inverted to give 

n 

kZ(i) = YK w fi® ■ (A27) 

13=1 



In terms of this new set of variables 

E (^«) 2 = E (e^(o) (e = E K«) 2 • (A28) 

We then have 

A = A(q) + A(s) + A(W L ) + A T (W) + A c (q, s, W L ) (A29) 

where 

A(q) = S Y,Vi 2 + 2 CWi • <?2 (A30) 
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A(s) = ]T si 2 + 2(V 4 C )si • s 2 



i=l 



(A31) 



A L (W) = E W L {i) 2 - 2(C£)# L (1) • W L {2) 



(A32) 



i=i 



where the longitudinal part of the tensor W is defined by 

n 

wT(0 = E*i*£(0 • 

The transverse contribution to A is given by 



(A33) 



/i=2 !/=l 



s (2) E(^;«) 2 - 2{C Q /R)w;{i)w;{2) 



i=i 



(A34) 



The term coupling the set <f, s, Wx, is given by 



A c (£ s, H> L ) = -2n5( 2 ) ]T £ • si + 2(gi • s 2 + g 2 • Si) V 2 C 



i=i 



-2(ps! + <%) • H> L (2) + 2(ps 2 + C' q 2 ) ■ W L (1) 



(A35) 



Notice that the transverse modes decouple from the longitudinal set coupled in A c . 

It will be very useful for us to rewrite A as a sum of a transverse part, already written 
down, and a logitudinal part which is a quadratic form in the vector 



K = q(l), q(2), s-(l), s(2), W L (1), W L (2) 



where we assume that the subscript a runs from 1 to 6. We have then 



(A36) 



Al — E 



(A37) 



where the matrix M is given explicitly by. 



M = 



So Co U\ u 2 W 4 

Co So u 2 ui — u 4 

ui u 2 C 4 u 3 

u 2 ui C 4 S (4) -u 3 

-« 4 -« 3 S( 2 ) C 2 

^W 4 M 3 C 2 
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(A38) 



Various quantities entering the matrix M are defined by 



a 



C 4 = V 4 C = j^Vtf(x) 



s (4) = v 4 q,U =0 = 8 



a 
1? 



C 2 = = -of 



Ml 



-nS® = - 



(A39) 
(A40) 
(A41) 
(A42) 



u 2 = V 2 C = aV 2 J(x) 



(A43) 



a 



«3 = -p=-^(V*/(x))' 



u 4 = -Cq = -Laf'(x) 



(A44) 



(A45) 



The complete change of variable from k to W in Eq. ( |A26|) requires noting that the 
Jacobian taking one from Ho IF is one. The argument of the exponential outside the 
average in Eq.( |A26| ) can also be written in terms of the set <p a and the transverse part of 
W . Then one has 



E 

i=l 



E W£« + * ■ S W 

/i,i/=i 



where 



6 2 n n 

i=l j/=1 /i=2 



(A46) 



h 1 = h 2 = 



(A47) 



/l 3 = 6(1) 

h* = 6(2) 
4 = 4(1) 



(A48) 
(A49) 
(A50) 
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he = t L (2) 



(A51) 



and 



(A52) 



We then have the result that G 2 factorizes into longitudinal and transverse components: 



G 2 ^,b) = G T (t T )G L (b,t L ) 



(A53) 



First consider the transverse contribution given by 
Gt^t) - 



2 n n dW"(i) 

mi n 



2vr 



cxp 



1=1 I/ = l |U=2 



j=l !/=l M =2 

This is a standard Gaussian integral which can be evaluated with the results 



n(n— 1) 



cxp 



EE 

M =2 i/=l 



EK(01 2 " 2/r^(l)C(2) 



where 



(A54) 



(A55) 



S T = S (2) 



(A56) 



Ct -~1i 



(A57) 



Ct 

St 



and 



(A58) 



7l = (1 - fi 



(A59) 



The longitudinal contribution to G 2 is also a standard gaussian integral: 



G L (b,t L ) = /n 



i=i 



iii (2tt)» 
1 



1 



(A60) 
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Thus the determination of G 2 reduces to an evaluation of the inverse and determinant of 
the matrix M given above. If we multiply M from left and right by the matrix 

1 \ 



Q = 



-1 I 
I 1 
-\ 1 
\ -1 

I 1 



(A61) 



then the new matrix 



M af3 = Y.QmQvP M i 



(A62) 



has the block diagonal form 



M 



V 



(So — C ) 





u x - u 2 





M 4 








\(So + Co) 





Ui + u 2 





■U4 


Ui - u 2 





\(S^-C 4 ) 





5U3 








Ux + u 2 





2(5^ + C 4 ) 





2m 3 


U4 










|(5 (2) +C 2 ) 
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2m 3 





2(5( 2 ) - 



(A63) 



It is easy to see that the inverse for M can be expressed in terms of the inverse of M as 

(M- 1 ) P = EW^(^" 1 W • ( A64 ) 

a/3 

Since M is block diagonal the evaluation of its determinant and inverse elements is 
straightforward. We find 



det M = D D E 

where Do is the determinant of the odd part of the matrix given by 
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(A65) 



2(S - C ) ui - u 2 u 4 

U 1 -U 2 \{S&-Ct) |«3 



«4 



|«s \{S^ + C 2 )j 



(A66) 



so 



2D Q = (So - Co) \(S^ - Ct)(S<» + C 2 ) 



u. 



— (u 1 — u 2 ) (ui - u 2 )(S ( - 2 ^ + C 2 ) - u 3 u 4 +u 4 u 3 (ui - u 2 ) - ■u 4 (S' (4) - C 4 ) . (A67) 



D E is the determinant of the even part of the matrix given by 



fi 



2 (S + C ) u Y + u 2 u 4 
Ul + u 2 2{S^ + C A ) 2u 3 
u A 2w 3 2{S^ - C 2 ) 



\ 



(A68) 



where 



D E /2 = (S + C ) \(S^ + C 4 )(^( 2 ) - C 2 ) 



u. 



- («i + u 2 ) f(«x + m 2 )(.S (2) - C 2 ) - M3M4I + m 4 [w 3 («i + "2) - u 4 (S {A) + C 4 )l (A69) 



Since it is easy to show that det Q = 1 we obtain 



det M = det Q det M det Q = det M = D E D Q 



(A70) 



The needed inverses are given by 



(M -1 )33 = 



(So-C )(S^+C 2 ) - 



Do 



(A71) 



/w-lN (^i - m 2 )m 4 - (So - C )m 3 
U w J35 — 



(A72) 



55 



(Sq-c )(s^ -c 4 )-( Ul -u 2 y 

Do 



(A73) 
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(M- 1 ) 



4d 



{S + C )(S^-C 2 )-uj 
D E 



(A74) 



(ui + u 2 )u A - (S + C )u 3 



(M ) 46 = 



(A75) 



(M _1 )66 = 



(g + Co)(,g( 4 ) + C 4 )-K+ M2 )- 



(A76) 



All odd-even inverse elements such as (M 1 ) 34 vanish and the rest of the elements follow 
using the symmetry of (M _1 ). One can then easily extract the inverse elements of (M _1 ) : 



(A77) 



(M" 1 ) 3 4 = (M- 1 ) 4 4- i (M- 1 )33 



(A78) 
(A79) 



(M-% = -(M- 1 ) 46 = \{M-%s - (M- 1 ) 46 



(M-^se = -(M" 1 ) 45 = -(M-^as + (M" 1 ^ 



(M- 1 ) 56 = ^(M- 1 ) 55 -(M- 1 ) 66 



(A80) 
(A81) 
(A82) 



In terms of / and its derivatives we find that we can write 



D F = 2a 3 



(1)12 



(A83) 



D = -a 3 
2 



(1)12 
O J 



^(M-^as = LV(1 - /) 



(0) 

o 



(A84) 
(A85) 



D (M- 1 )35 = L C r 2 (l-/) K 



(i) 
o 



(A86) 
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Do(M-% b = -a 2 4 ] (A87) 
D E (M-% A = L 2 a 2 (l + /)4 0) (A88) 



D E {M~ 1 ) m = La 2 (l + (A89) 

D E (M~ l ) m = a 2 4 ] . (A90) 

These results give the explicit L dependences of the various matrix elements. The k's are 
independent of L and given by 

«g> = (V 2 /)' - (A91) 

4 2) = (l + /)(V 4 /) + -[(V 2 /) + ] 2 (A92) 

4 0) = (/")- " (A93) 

"& = (V 2 /)' + ^^^^ (A94) 

4 2) = (l-/)(V 4 /)- + [(V 2 /)_] 2 (A95) 

4 0) = -(/")+ " YZJ (A96) 
where we have introduced the notation 

A± = A(x) ± A(0) . (A97) 
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TABLES 



„(1) 

K E 


2x 





K E 


16 


2n 


(0) 
K E 

K (0) 

Kq 

K (l) 
Kq 

K (2) 
Kq 


x 2 

12 

X 3 

2 

-2x 4 


1 


-2n 
1 


D e 


16a 3 x 2 


4racr 3 


Do 


a 48 


n<r 3 


del M 


3 


4nV 6 



TABLE I. small x (left) and large x values for various quantities (far left) denned in the text 
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FIGURES 




12 3 4 

X 



FIG. 1. The most probable scaled velocity of vortex 2 multiplied by the magnitude of the 
scaled distance of separation x, u = x ■ u{2)L/Tc versus x. This velocity is directed along x the line 
connecting the two tagged vortices. The most probable velocity of vortex 1 is equal and opposite 
to that of vortex 2. 
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u 



FIG. 2. Probability distribution (unnormalized) for the scaled longitudinal velocity 
u = x ■ u(2)L/Tc = —x ■ u(l)L/Tc in the small scaled distance x limit. The peak in this curve gives 
the small x limit of the quantity plotted in Fig.l. 
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-2 2 4 6 8 10 

U 

FIG. 3. Plot of unnormalized probability P[vi,V2] for different values of the scaled dis- 
tance x between the two tagged vortices versus the scaled velocity in the center of mass 
u = x ■ u(2)L/Tc = —x ■ u(l)L/Tc. The normalization changes with x so the different heights 
of the curves is not significant in this plot. The curves, as one moves from left to right, are labelled 
by x = 3.0, 2.0, 1.0 and 0.7 respectively. 



50 




51 



